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Isotropic soft-core potentials with two characteristic length scales have been used since 40 years 
to describe systems with polymorphism. In the recent years intense research is showing that these 
potentials also display polyamorphism and several anomalies, including structural, diffusion and 
density anomaly. These anomalies occur in a hierarchy that resembles the anomalies of water. 
However, the absence of directional bonding in these isotropic potentials makes them different 



from water. Other systems, such as colloidal suspensions, protein solutions or liquid metals, can 
be well described by these family of potentials, opening the possibility of studying the mechanism 



generating the polyamorphism and anomalies in these complex liquids. 
PACS numbers: 64.70.Ja, 82.70.Dd, 61.20. Ja, 64.70.qj, 65.20.De 



I. INTRODUCTION 

Isotropic pair interaction potentials are usually considered the prototype for simple atomic 
systems, such as argon. The most famous among them is the phenomenological potential 
proposed in 1931 by Sir John Edward Lennard- Jones (LJ)i, commonly adopted as the text- 
book model for real gases. The LJ potential incorporates the short-range repulsion, due to 
the Pauli's quantum exclusion principle among electron orbitals, as a function ~ 1/r 12 of 
the distance r between the centres of mass of the atoms. It also includes the van der Waals 
attraction, due to instantaneous induced dipole-dipole London dispersion forces between the 
electron clouds, as a long-range function ~ —1/r 6 . These two components are enough to 
generate a phase diagram with a gas, a liquid and a solid phase, as for neutral atoms or 
simple molecules. The LJ model reproduces not only the thermodynamics, but also the 
dynamics and the kinetics of these systems, providing for example a good starting point for 
studying processes such as the homogeneous nucleation of the crystal phase. 

The L J model, and similar potentials such as square wells, are useful also for more complex 
systems, e. g. colloidal suspensions or protein solutions. In these cases these isotropic pair 
interaction potentials can be used to represent the interactions between the particles of the 
solute when the degrees of freedom of the solvent are implicitly taken into account in the 
effective interaction potential. However, there are (anomalous) properties of these and other 
systems, e. g. liquid metals or water, that cannot be reproduced by simple potentials. It is, 
therefore, natural to ask if the family of isotropic potentials can be extended in such a way 
to describe the phase diagram of the anomalous substances. 

Anomalous systems such as water and silica are network-forming liquids with strongly 
anisotropic interactions. However, for other systems such as liquid metals 2-14 , colloids^— or 
biological solutions^ - — the use of soft-core isotropic potentials with two characteristic length 
scales is a particularly suitable way of constructing effective pair interactions capable of 
describing the anomalies of these systems. Without the pretension to complete or exhaustive 
coverage, here we recall some recent results about this topic. For other aspects related to 
soft-core potentials we remit to previous reviews 2 ^"—. 



II. ANOMALOUS LIQUIDS 

Experiments for Ga 7 , B&, Tb21, S 2 ^, Be, Mg, Ca, Sr, Ba 30 , Si0 2 , P, Se, Ce, Cs, Rb, 
Co, Gei2, Gei 5 Te 8 g^i and simulations for SiOf^— , S^ and BeF^ reveal the presence of 
a temperature of maximum density (TMD) at constant pressure below which the density 
decreases when the temperature is lowered isobarically. This behavior is at variance with 
that of normal (or argon-like) liquids where the density monotonically increases when the 
temperature is decreased at constant pressure. The most famous example of liquid with this 
anomaly is water, whose TMD at 1 atm is at 4°C. Below the TMD the isobaric thermal 
expansion coefficient ap = y f^L of water assumes negative values. In normal liquids ap is 
always positive because it is proportional to the (always positive) cross-fluctuation of volume 
and entropy. Other thermodynamic anomalies of water include the anomalous increase 
of isothermal volume fluctuations, proportional to the isothermal compressibility Kf = 
y §75 L, below 46°C at 1 atm and the anomalous increase of isobaric entropy fluctuations, 
proportional to the isobaric heat capacity Cp = T |^| , below 35°C at 1 atm 37 . 

Another anomaly observed in water and other liquid is related to the diffusion coefficient 

D, defined as 

D , lim { -^fl (!) 
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where t is the time, d = 3 is the dimension of the system, 

(Ar(t) 2 )^([r(t + t)-r(t )] 2 ) (2) 

is the mean square displacement of a single particle, t is any time at equilibrium and 
the average is over the initial t and over the particles in the system. In a normal liquid 
D decreases when density p or pressure P are increased. Anomalous liquids, instead, are 
characterized by a region of the phase diagram where D increases when P is increased at 
constant temperature T. In the case of water, for example, experiments show that the 
normal behaviour of D is restored only at pressures higher than P « 1.1 kbar at 283 K— . 

Also the structure can be anomalous. Normal liquids tend to become more structured 
when compressed. This can be quantified by a translational order parameter t that measures 
the tendency of the molecules to adopt preferential separations, and by an orientational order 
parameter Qi that measures the tendency of a molecule and its nearest neighbours to assume 
a specific local arrangement, as considered by Steinhardt et al.— . The translational order 



parameter is defined as 3 ^ 4 ^ 41 - 
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where £ = rp 1//3 is a reduced distance (in units of the mean interparticle separation p -1 ^ 3 ) and 
g(£) is the radial distribution function. As the parameter t depends only on the deviations 
of g(£) from unity, its value is sensible to long range periodicities. For an ideal gas g(£) 
is constant and equal to 1 and there is no translational order (t — 0). For a crystal phase 
g(£) 7^ 1 for long distances and t becomes large. 

The orientational order parameter is by definition 39 
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where 1 = 1,2, ... is an index, Q\ = ^ YZ=-i \( Y L)k 

neighbour particles and (l^Jfc — \ Ej=i^m( r !i) * s the average of the spherical harmonics 
Y\ m with indices I and m, evaluated over the vector distance ry between particles z and j. 
For / = 6 and fc = 12, Q6 reaches its minimum value Qg h = \j\fk = 0.287 for an isotropic 
homogeneous system, while for a fully ordered f.c.c arrangement is Qg cc = 0.574. 

For a normal liquid t and Qz increase with pressure. Anomalous liquids, instead, show a 
region where the structural order parameters decrease for increasing pressure (or density) at 
constant T, i.e. the system becomes more disordered. This is what has been observed, for 
example, in molecular dynamics simulations for water by Errington and Debendetti^ and 
by Shell et al. for silica 34 . 

All the anomalies of water have a well determined sequence found in experiments 3 ^ and 
simulations 4 ^. In the T — p plane the water structural anomaly region is encompassing 
the diffusion anomaly region which includes the density anomaly region. However in other 
liquids the sequence of anomalies may be different. For example, for silica the anomalous 
diffusion region contains the structural anomaly regions that, in turn, includes the density 
anomaly region^ 4 -. 

A. Polymorphism, Polyamorphism and Liquid liquid phase transition 

Another anomaly that has received considerable attention in recent years is the possible 
existence of a liquid-liquid (LL) phase transition for single- component systems with a stan- 
dard gas-liquid critical point. The two coexisting liquids, the high density liquid (HDL) and 



the low density liquid (LDL), would differ in density and local structure, as proposed by 
Poole et al. in 1992 42 , based on simulations for a water model. At that time it was known 
that water can have more than an amorphous stated, a high density amorphous (HDA) 
and a low density amorphous (LDA), separated by discontinuous density-change with the 
character of first-order phase transitions. An even higher amorphous state, the very HDA 
(VHDA), of water has been observed in 2001 by Loerting and coworkers 44 . 

This property, called polyamorphism, i.e. the occurrence of more than one amorphous 
state, is usually associated to the polymorphism, i.e. the occurrence of more than one crystal 
phase 45 . Typical examples of polymorphic substances are water with at least 16 forms of 
ice, and carbon with diamond and graphite (made of graphene sheets, as are fullerenes and 
carbon nanotubes). 

Direct evidences for liquid polyamorphism have been observed experimentally in 
phosphorous 4 ^— in 2000, triphenyl phosphite 4 ^ - — in 2004 and in yttrium oxide-aluminum 
oxide melts^ in 2008. Experimental data consistent with a LL phase transition have also 
been presented for single- component systems, besides water, such as silica^ 5 - 15 "^ 5 - 4 -, carbon 5 - 5 -, 
selenium 5 - 6 -, and cobalt 57 , among others 5 -^— . 

A LL critical point has been predicted by simulations for all the commonly used mod- 
els of water, including SPC/E, ST2, TIP4P and TIP5P— >§i, and for specific models of 
phosphorous^, supercooled silica^ 5 -^' 5 ^^, hydrogen^ 4 -, and carbon based on molecular dy- 
namics simulations 6 - 5 -, For carbon, however, subsequent ab initio simulations 66 and simula- 
tions for a semiempirical potential partly based on ab initio data 6 - 7 - did not confirm this 
finding. 

III. ISOTROPIC MODELS WITH ANOMALIES 

The anomalies described in the previous section can be reproduced by using isotropic 
core-softened potentials. The core-softened potentials are usually characterized by a change 
of curvature within the repulsive range, such as a ramp or a shoulder. 

• Ramp-like potentials have their softened region defined by a repulsive ramp, establish- 
ing two competing equilibrium distances, and in some cases an attractive well 6 ^"— . 

• Shouldered potentials are composed by a hard-core, the repulsive shoulder softening 



the core, and an optional attractive well 1 



Thermodynamic anomalies have been widely reported for such potentials. Some of these 
anomalies depend on the details of the potentials, as we will show in the following. Among 
these potentials, those with an attractive part display more than one first-order phase tran- 
sition. 

It was after almost 40 years from the introduction of the LJ potential that Hemmer and 
Stell explored the idea of core-softening, often referred as core-collapse, proposing the first 
repulsive-ramp potential^ 1 ^. They realised that by making the core more penetrable they 
could induce a second first-order phase transition in systems that already have one. They 
justified the occurrence of the second transition in a lattice gas by means of a heuristic 
argument based on the (particle-hole) symmetry between occupied and unoccupied cells. 
This argument, however, is not applicable to continuous systems where the particle- hole 
symmetry does not hold. Nevertheless, they showed explicitly for fluids in one dimension 
how a similar phenomenon could appear by softening the hard-core, finding a range of 
parameters for which the phase diagram displays two first order phase transitions, both 
ending in critical points. The kind of potentials they proposed were discontinuous. However, 
they argued that other analytical potentials arbitrarily close to those they considered could 
be constructed and the second transition would persist for such potentials. The high-density 
critical point was first interpreted as a solid-solid isostructural phase transition for systems 
such as Cs or Ce^2. In their original work— ^ Hemmer and Stell remarked that for one 
dimensional models with long range attraction the isobaric thermal expansion coefficient otp 
can take anomalous negative values. 

In 1976, Kincaid et al.— && proposed a variant of the Hemmer-Stell potential with the 
hard core softened by adding a finite shoulder of constant positive energy. They found 
at high density and pressure an isostructural phase transition. Kincaid and Stell^ extended 
this work by introducing a shouldered potential without attraction to describe solid mixtures 
with isostructural transitions. 

After this pioneering works, several soft-core potentials have been proposed and analyzed 
with approximate methods or numerical simulations to study the properties of complex 
liquids as liquid metals, alloys, electrolytes, colloids and, to some extent, water-"— liLikiSiriSi 



For example, Stillinger and Weber— in 1978 used molecular dynamics simulations to 
study the phase behaviour of the Gaussian core model. They reported the surprising and 
unexpected finding that the model at equilibrium displays water-like anomalies, such as a 
density of maximum melting temperature, hence a region of decreasing volume upon melting, 
a negative thermal expansion coefficient in the fluid phase and an increase of self-diffusion 
upon isothermal compression. 

Later, in 1991, Debenedetti et al.— showed that anomalies can occur in two-length-scales 
potentials also when the inner distance is attractive and the larger distance is repulsive. 
They showed that such a potential on a lattice can induce the formation of an open (low 
density) structure at low P and T. Upon heating or pressurization the open structure looses 
stability and collapses into a closed (denser) structure, leading to the anomaly in density 
and to negative op. A similar model on a lattice was studied by de Oliveira and Barbosa^ 
finding a LL coexistence with a line of critical points. More recently Archer and Wilding 106 
studied a different isotropic model with short range attraction and larger range repulsion, 
supporting the existence of a similar line of critical points. 

In 1993 Head-Gordon and Stillinger— »^ showed that from the inversion of the radial 
distribution function of water is possible to deduce an effective oxygen-oxygen interaction 
potential that resembles a continuous version of the Stell-Hemmer shouldered potential. 
In 1996, Cho et al . 107 ' 108 extended even further the original idea of Stell and Hemmer of a 
potential with two length scales, by proposing a potential with an inner well and an outer well 
with two characteristic attractive energies separated by a local maximum. They analyzed 
the model in one dimension showing the presence of the density anomaly and proposed to 
consider it as an effective potential for the second shell of water, assuming that the first 
shell can be considered as part of an invariant inner core that does not play a relevant role 
in the anomalous density behavior of water. 

In 1998, Sadr-Lahijany et al.— ^ considered a similar shouldered potential in two dimen- 
sions, both in the discrete and the continuous version, with a deep attractive well. The 
shoulder gives rise to an inner distance that is less attractive than the outer distance. They 
found polymorphism: a low-P triangular lattice less dense than the liquid, and a high-P 
square lattice denser than the liquid. Furthermore, they reported the occurrence of three 
anomalies: the density anomaly, the increase of isothermal compressibility upon cooling, and 
the diffusion anomaly. By a simplified argument in one dimension, it was shown how the in- 



crease of P reduces the Gibbs free energy at the inner distance of the potential, with respect 
to the larger and more attractive distance, inducing the collapse to the denser structure. 
The interplay between the two different local structures near the freezing line was proposed 
as the mechanism responsible of the anomalies. This rationalization was better clarified by 
Wilding and Mageo^. Further studies on a continuous version of the model by Netz et al.— 
showed anomalous behavior in the stable region of the phase diagram if the outer minimum 
is deeper than the inner minimum. In the case of a deeper inner minimum, anomalous 
behavior occurs inside the unstable region. Nevertheless, the study of the same model in 
three dimensions at low T in a stable liquid state by Quigley and Proberti 10 suggests that 
the anomalous behavior of such model is unique to the two-dimensional case. 

Also in the 1998, Jagla 70,71 studied a version of the purely repulsive Kincaid and Stell 
potential by softening the core with a linear ramp, as in the original Hemmer and Stell 
potential. The absence of the attractive term implies the absence of the liquid-gas phase 
transition. Under pressurization the potential displays many crystalline polymorphs and 
anomaly in density. Also in this case the occurrence of these anomalies is associated to the 
competition between the two scales of the potential: the hard-core distance and the distance 
at which the interaction becomes zero. The appearance of complex ground state structures 
in the system is a consequence of the competition between two terms in the enthalpy H: the 
pressure-volume term PV leads to minimize the volume, while the soft-core repulsive term 
leads to maximize the interparticle distance. The different arrangements of particles in the 
ground state depend on the values of the pressure P and the values of the two length scales 
of the potential. This result was later confirmed by Kumar et al.— using integral-theories. 
Polymorphism for purely repulsive core-softened potentials was also obtained a few years 
later by Velasco et al. in collaboration with Hemmer and Stel l 141111 . More recently, by using 
free-energy calculations, Gribova et al.— showed that for this potential the system exhibits 
the water-like density and diffusion anomaly and that the anomalies move to the region 
where the crystal is stable with increasing repulsive-step width. 

In 1999, Jagla analyzed a slightly modified version of the Hemmer and Stell potential 7 - 2 - 
in three dimensions showing by simulations the appearance of water-like anomalies and 
predicting by theory a LL critical point. However, in the simulations the LL critical point 
results inaccessible due to the inevitable crystallization 72 . 



In 2001, by adopting different parameters for the potential^ and by analyzing a con- 
tinuous version of the model 74 , Jagla showed that the LL phase transition is observed in 
simulations in two and three dimensions, when the attraction is strong enough, while dis- 
appears if the attraction strength is approaching zero^ 4 -. In both case, however, the density 
displays hysteresis at low T due to mechanical metastabilities, and water-like anomalies are 
observed, regardless the absence of the LL critical point, showing that the presence of the 
anomalies is not a sufficient condition for the occurrence of the LL critical point^ 4 -. Later, 
Wilding and Magee observed that the locus of density maxima ends very close to the LL 
critical point, when this is present^. Xu et al . 112 ' 113 showed that in the super-critical region 
of the LL critical point the dynamics changes from Arrhenius to non-Arrhenius when the 
line of maxima correlation length, approximated with the line of maxima of the specific 
heat Cp, is crossed. Caballero and Puertas^ 4 - studied the model by first-order perturbation 
theory for different choices of the attractive range, finding that the LL phase transition is 
attraction-driven for long ranged potentials, and is compression-driven when the interaction 
is shortened. 

In the same year 2001, Franzese et al.— showed that the occurrence of the anomalies 
is not a necessary condition for the LL critical point. They found, by simulations and 
integral equation theory^, that a discontinuous shouldered well (DSW) potential in three 
dimensions, similar to that originally proposed by Kincaid et al., displays a LL coexistence, 
ending in a critical point, metastable with respect to the crystal. They also showed that 
the system does not display an anomalous behavior in density. This result suggested that 
the type of systems displaying the LL phase transition could be broader than what was 
previously hypothesized, and that experimental evidences of a LL phase transition should 
be seeked also in systems without density anomaly. In 2004 Skybinsky et al.— showed that, 
by changing the parameters of the DSW, the phase diagram displays a LL critical point 
that is stable with respect the spontaneous homogeneous crystallization. They developed a 
modified van der Waals equation that qualitatively reproduces the behavior of both liquid- 
gas and LL critical points of the model. In 2005 Malescio et al. extended the previous 
work to the case of large soft-core ranges, by using an integral equation approach in the 
hypernetted-chain approximation. They showed that only a limited range of parameters of 
the DSW model give rise to a phase diagram with an accessible LL critical point and that this 
occurs when the repulsive component of the potential equilibrates the attractive component, 



in particular when the repulsive volume weighted by the repulsive energy compensates the 
attractive volume weighted by the attractive energy 9 -^. Lattice Monte Carlo simulations 
by Balladares and Barbosa 9 ^ for the DSW confirmed the LL coexistence but gave a phase 
diagram with a line of critical points connecting the LL phase transition and the liquid- 
gas phase transition. This feature appears to be an artifact of the lattice, because off- 
lattice Monte Carlo simulations by Rzysko et al.— showed no such a line of critical points, 
confirming the previous results of molecular dynamics simulations and theory 9 ^ ^ 115 ' 116 and 
providing an accurate estimates of the LL critical points and its exponents. 

In 2003, Buldyrev and Stanley 92 tested the idea that for any characteristic length scale 
in the interaction potential the system possibly displays a different liquid phase and a new 
LL phase transition. They added an extra discontinuous step to the soft-core of the DSW, 
defining a third characteristic distance. For certain values of the parameters of the potential 
the system presents up to three first-order phase transitions between fluids of different 
densities ending in critical points. The radial distribution functions g(r) displays dramatic 
differences in structure between the low density liquid (LDL) and the high density liquid 
(HDL), but less pronounced differences between the HDL and the very high density liquid 
(VHDL). These results suggest that more critical points could be created by adding more 
steps to the potential and carefully selecting the parameters. However, for ksT larger than 
the steps of the potential, the effect of the new length scales becomes negligible and the 
phase diagram converges to that of a system with a continuous potential. 

Nevertheless, a few years later Cervantes et al.— showed, by calculating the free energy 
by discrete perturbation theory, that there is a range of combination of parameters of the 
two-scales DSW that generate a phase diagram with three critical points: between gas and 
liquid, between LDL and HDL, and between HDL and VHDL. The three critical points for 
the two-scales DSW were found also, by Artemenko et al.— , by analyzing the modified van 
der Waals equation proposed in Ref.— . Therefore is not necessary to add more than two 
length scales in the potential to reproduce phase diagrams with more than two fluid-fluid 
critical points. It is enough to have two characteristic length scales that compete creating 
a multiplicity of minima for the free energy at different P and T. These results, about the 
minimum conditions to reproduce polyamorphism including at least four different charac- 
teristic densities, are particularly interesting considering that experiments^ and molecular 
models simulations^ show the existence of a very high density amorphous phase of water—. 

10 



Polyamorphism and polymorphism were also studied in two dimensions in the Kincaid- 
Stell purely repulsive shouldered model by Malescio and Pellicane in 2003^ and in a purely 
repulsive potential with a shoulder and a long-range repulsive tail by Camp^^. The fluid 
phase at low temperature exhibits a very rich variety of structures, including chains, stripes 
and polygons 122 . The Camp model approximates a two-dimensional system of dipolar par- 
ticles in a strong field aligned perpendicular to the pian o 123 ' 124 . 

The properties of the core-softened potentials depend on the ratio between the two charac- 
teristic scales. Yan and co-workers 7 ^* 7 ^ explored a range of systems going from a hard-sphere 
potential to a pure ramp without hard-core, finding that thermodynamic (negative op) and 
dynamic (diffusion) anomalies occur almost across the entire range, while water-like struc- 
tural anomalies occur only for cases with the ratio between the two length scales comparable 
with that between the first two peaks of water in standard conditions. They also showed 
that the anomalies have the same hierarchy as in water: the density anomaly occurs within 
the region of diffusion anomaly, that is found within the region of structural anomaly, as 
observed for the SPC/E water—. This similarity with water exists despite the lack of di- 
rectionality in the isotropic potentials. The analogy with water was pushed forward by 
Yan et al. 79 comparing the anomalies in the super-critical region of the two scale potential 
with those of TIP5P water, giving new support to the initial proposal by Cho et al.— that 
two-scale isotropic potentials could represent an effective interaction for the second shell of 
water, being the first part of the core. 

The water-like hierarchy was found also by de Oliveira et al.— ^ in three dimensions for 
a shouldered potential with a small attractive region, that resembles the Camp potential 95 
and that has been extensively studied by Pizio et al. recently 125 . These results on all these 
soft-core models show that orientational interactions, such as hydrogen bonding, are not 
a necessary condition for the presence of water-like anomalies. The anomalous features at 
about the temperature of the maximum density are caused by the reduction of the large 
empty spaces around the molecules upon compression and heating. This effect is captured 
by the soft-core potentials with two characteristic length scales. However, other anomalies 
of water at lower T, such as the liquid-liquid phase transition and the associated Widom line 
are not reproduced in the correct way. In particular their slope in the P — T phase diagram 
is opposite to what is observed in molecular water models such as TIP4P. This difference 
implies an opposite behaviour as function of P for thermodynamic and dynamic quantities 
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in the supercooled regime. 

All the example of potentials mentioned above have been proposed for systems going from 
water to colloids, to protein solutions, to liquid metals, with specific experimental cases, 
such as polystyrene monolayers between water and air in two dimensions^ 7 -, or gallium 7 -, or 
tellurium^ or liquid alkaline-earth metals near the melting points in three dimensions. 
The effective (pseudo)potentials representing these systems are often similar in shape, but 
with different details. It is therefore interesting to understand how their properties depend 
on the detailed features of the potential. To this goal we present in the next section a case 
study for a highly tunable potential that offers the opportunity to explore the properties of 
a large family of soft-core potentials. This potential is a continuous version of the DSW, 
the continuous shouldered well (CSW) potential, introduced by Franzese in 2007 97 . It was 
shown that the CSW potential has density anomaly in three dimensions 97 and a detailed 
analysis by de Oliveira et al.— revealed also diffusion anomaly and structural anomaly with 
a water-like hierarchy. Recently Standaert et al.— have adopted this model to study the 
condition of anomalous (non- Gaussian) self-diffusion in a system driven out of equilibrium 
by intermittent length rescaling. 

IV. THE CONTINUOUS SHOULDERED WELL POTENTIAL 

The CSW model^ 7 - consists of a set of identical particles interacting through the isotropic 
pairwise potential 



U(r) = 777- . . . — Ua exp 

KJ l + exp{A{r-R R )/a) 
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(Fig. [T]) where a is the diameter of the particles, Ra and Rr are the distance of the attractive 
minimum and the repulsive radius, respectively, Ua and Ur are the energies of the attractive 
well and the repulsive shoulder, respectively, 5\ is the variance of the Gaussian centered in 
Ra, and A is the parameter which controls the slope between the shoulder and the well at 
Rr. Varying the parameters the potential can be tuned from a repulsive shoulder to a deep 
double well. In particular, by increasing A the soft-core repulsion becomes more penetrable 
near the minimum of the attractive well, and the softness of the potential increases for 
r > R R and decreases for r < R R . We fixed here the set of values Ur/Ua = 2, -Rr/<2 = 1.6, 
Ra/o, = 2, {5a/ o) = 0.1, while we change A, considering the values A = 15, 30, 100, 300, 500 
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FIG. 1: The Continuous Shouldered Well (CSW) potential for A = 15, 30, 100, 300, 500 (continuous 
lines) and the Discontinuous Shouldered Well (DSW) potential (dotted black line). By increasing 
A the CSW potential approximates the DSW around Rr. 

(Fig. [1]) going from the case A = 15 studied in Ref.— ^ to slopes that approach the infinite 
value of DSW. 

V. THE PHASE DIAGRAM 

The phase diagram for all the considered values of A display the same qualitative be- 
haviour in the plane P* — p* (Fig. |2]) where P* = Pa 3 /e and p* = pa 3 are reduced pressure 
and density, or P* — T* (Fig. [3]) where T* = ksT/e, with ks Boltzmann constant, is the 
reduced temperature. In particular, at low T* the isotherms are non-monotonic (van der 
Waals loops), corresponding to the coexistence of (i) gas and liquid at low p*, (ii) low density 
liquid (LDL) and high density liquid (HDL) at higher p*. The two coexistence regions end 
in critical points: C\ for gas-liquid coexistence, C2 for LDL-HDL coexistence. For A < 100 
the HDL phase is metastable with respect to the crystal, but with a lifetime long enough to 
allow us to equilibrate the liquid around CV For A > 100 the HDL lifetime is too short to 
equilibrate the liquid around C 2 and we extrapolate the location of C 2 by linear fitting from 
higher T* isotherms. 
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FIG. 2: Isotherms in the P* — p* phase diagram from simulations for CSW potential for A = 30 
(a), A = 100 (b), A = 300 (c) and A = 500 (d) (in all the panels, from top to bottom, T* = 1.4, 
1.3, 1.2, 1.1, 1.0, 0.9, 0.8, 0.7, 0.6, 0.55, 0.5 and 0.45). 

As A is increased, both C\ and C<i tend to the corresponding values for the DSW potential. 
This result is consistent with the idea that the DSW can be seen as a limiting case of the 
family of CSW potentials presented here. However, temperature and pressure of the LDL- 
HDL critical point for A = 500 (T£ 2 = 0.52 ±0.01, P£ 2 = 0.204 ±0.007, p* C2 = 0.272 ±0.008) 
are still far from the values for the DSW potential (T£ 2 = 0.69 ± 0.02, P£ 2 = 0.110 ± 0.002, 
p* C2 =0.280 ±0.020). 

VI. ANOMALIES 

The CSW has anomalies in density, diffusion and structure. By numerical simulations, de 
Oliveira et al.— studied the case A = 15 finding the same hierarchy of anomalies as reported 
for other two-scale potentials and for the SPC/E water. 

In Ref.— we analyze in details the cases A = 30, 100, 300, 500. In summary, for all the 
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FIG. 3: The P* -T* phase diagram for A = 30 (a), A = 100 (b), A = 300 (c) and A = 500 (d). In 
all the panels, we show (from bottom at T* = 1.2) isochores from p* = 0.04 to 0.20 (with constant 
density separation 5p* = 0.01), from 0.205 to 0.215 (with 5p* = 0.005), and from 0.22 to 0.25 (with 
5p* = 0.01). The points where isochores cross correspond to the gas-liquid and LDL-HDL critical 
points {C\ and C*2, respectivelly). The black dashed line at low T is a guide for the eye estimating 
the TMD line, correspondig to the line of minima along the isochores. 

considered values of A, we observe the following. 

• Density anomaly. 

We find a temperature of minimum pressure (dP/dT) = along the isochores near 
the LDL-HDL critical point Ci- These temperatures correspond to the TMD line 
at constant P, because the condition (dP/dT) v = implies, according to Maxwell 
relations, (dS/dV) T = 0. Consequently, it is (dS/dP) T = (dS/dV) T (dV/dP) T = 
and, again using Maxwell relations, (dV/dT) p = (dp/dT) p = 0, which implies the 
presence of the temperature of maximum density (TMD) at constant P. 
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Diffusion anomaly. 

We find an anomalous-diffusion region, i.e. a region (pDmin < P < PDmax) where D, 
defined by Eq.flT}, increases with increasing density at constant T. 

Structural anomaly. 

We find that t, defined as in Eq.([n]) increases with increasing density, for p < pt max , 
and reaches a maximum at pt m ^- Above p< max , for increasing p, t decreases until it 
reaches a minimum at pt . . For p > p t . , t recovers the normal behaviour. 

Moreover, we observe that Qq, as defined in Eq.fi]), has a non-monotonic behaviour 
along the isotherms with a maximum at PQ max . The density pQ m ^ for each isotherm lies 
between pt max and pt min . In the area between p<g max and p tmin both order parameters 
decrease for increasing p, hence the liquid becomes more disordered for increasing 
density. This behaviour defines the structural anomaly region (pQ max < p < Pt min )- 

Hierarchy of anomalies. 

By varying the values of A, the regions of anomalies are affected in different ways 
while the hierarchy in which they occur is always preserved (Fig. H]). We find that 
the anomalous region where p decreases for decreasing T shrinks for increasing A and 
possibly tends to collapse onto one single point in the P* — T* plane for A — > oo. This 
result would be consistent with the behaviour of the DSW potential, that does not 
show density anomaly^" ^ 1151116 . Also the region of the diffusion anomaly contracts as 
A increases. We find that the diffusion anomaly region always encompasses the TMD 
line between p Dmin and p Dm ax (Fig. SJ). 

The region of structural anomaly does not contract by increasing the value of A 
but tends asymptotically to a fixed region in the T* — p* plane (Fig. 2]). This weak 
dependence on A suggests that the occurrence of the structural anomaly does not 
disappear for very steep soft-core potentials. This prediction is consistent with the 
excess entropy calculations that allow de Oliveira et al.— to argue that the structural 
anomaly should be observable also for the DSW potential, here considered as the limit 
of the CSW for A ^ oo. 
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0.25 # 0.2 
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FIG. 4: Hierarchy of anomalies in the T* — p* plane, plotted for (a) A = 15, (b) A = 30, (c) 
A = 100 and (d) A = 500. As the value of A increases, both the regions of density anomaly and 
diffusion anomaly contract, with the diffusion anomaly region always encompasses the TMD line. 
The region of structural anomaly is only weakly affected. 

VII. OUTLOOK 

Isotropic soft-core potentials with two characteristic distances are able to display the 
complex behaviour of anomalous liquids. They can describe the effective interactions of 
systems such as colloids, protein solutions or liquid metals. Due to the lack of directional 
interactions, these potentials provide a mechanism for the anomalies that is alternative to 
the bonding of network-forming liquid, e. g. water. As a consequence, their use as coarse- 
grained models of water is seriously disputed. In particular, they have a P-dependence of 
the structural fluctuations and a supercooled phase diagram that are different from those 
found in simulations of molecular models of water, such as SPC/E, ST2, TIP4P, TIP5P 
(see for example^) or of other coarse-grained water models with explicit hydrogen bond 
interactions*^-—. This wrong P-dependence affects, not only the supercooled state, but 
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also properties at ambient T, such as the velocity of sound. Water sound propagation 
shows a discontinuity in P that has a negative slope in the P-T phase diagram at 293 K 
and 0.29 GPa^, consistent with the slope of the maxima of structural fluctuations found 
at the same T and P in numerical simulations of TIP4P-water by Saitta and Datchi^. 
Nevertheless, the isotropic soft-core potentials with two characteristic distances display a 
sequence of anomalies that is water-like. 

To elucidate their properties we have presented here the case of the CSW potential 
whose anomalous behaviour is affected by the parameter A associated to the steepness of 
the repulsive shoulder. For all the considered values of A, the phase diagram displays two 
first-order phase transitions, corresponding to a liquid-gas phase transition at low densities 
and a liquid-liquid phase transition at higher densities, both ending in critical points. For 
15 < A < 500 we verify that the anomalies in density, diffusion and structure are in the 
same hierarchy as in water. We show that, as the value of A increases, the regions of density 
and diffusion anomaly contract in the T — p plane, while the region of structural anomaly 
is weakly affected. 

Since by increasing A the CSW potential approaches the discontinuous shouldered well 
potential (DSW), the contraction of the density anomaly and diffusion anomaly for A — > oo 
is consistent with the fact that for the DSW potential no TMD is observed. Our results 
suggest that the structural anomalies, instead, should be present also in the DSW poten- 
tial. Therefore, the isotropic soft-core potential with two characteristic distances can be 
considered as the prototype for anomalous liquids with no directional bonding. The fact 
that its properties depend in a dramatic way on the details of the potential calls for further 
investigations to understand better the mechanisms that regulates the appearance of the 
anomalies. 
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